Interaction between Kondo impurities in a quantum corral 
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We calculate the spectral densities for two impurities inside an elliptical quantum corral, using 
exact diagonalization in the relevant Hilbert subspace and embedding into the rest of the sys- 
tem. Fore one impurity, the space and energy dependence of the change in differential conductance 
AdI jdV observed in the quantum mirage experiment is reproduced. In presence of another impu- 
rity, AdI jdV is very sensitive to the hybridization between impurity and bulk. The impurities are 
correlated ferromagnetically between them. A hopping > 0.15 eV between impurities destroys the 
Kondo resonance. 
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In recent years, the manipulation of single atoms on 
top of a surface using scanning tunneling microscopy 
(STM) was made possible, and quantum corrals have 
been assembled by depositing a closed line of atoms or 
molecules on noble metal surfaces . [^]-[| The local con- 
duction spectral density of states p c (r,u>), measured by 
differential conductance dl/dV reveals patterns that re- 
mind the wave functions of two-dimensional noninteract- 
ing electrons under the corresponding confinement po- 
tential. In a recent experiment, a Co atom has been 
placed at a focus of an elliptic quantum corral, and the 
corresponding Kondo feature is observed not only at that 
position, but also at the other focus, where a "mirage" is 
formed as a consequence of the quantum interference. 
Several variants of this experiment, some of them involv- 
ing several impurities (Co atoms) and eventually mirages 
inside the corral, are being performed. Q The main fea- 
tures of the observed space and voltage dependence of 
dl/dV have been reproduced by several theories. H 
Since in Refs. the density of states per spin at the 

impurity Pdiy>) is assumed rather than calculated, these 
theories cannot account for the interaction between im- 
purities. In Ref. H the Kondo effect is absent, and per- 
turbation theory in the Coulomb repulsion U 0,0 is 
restricted to small values of U . 

The aim of the present work is to present a theory 
of the quantum mirage which is able to reproduce the 
experimental results for the case of one impurity and 
give reliable predictions when more than one impurity 
is inside the corral. We show that experiments with two 
impurities can elucidate the role of the direct hybridiza- 
tion between the impurity and the bulk Vb- Scattering 
theories obtained an excellent agreement with ex- 

periment assuming that the resonant level width due to 
hybridization with bulk states 6b is as large as that due 
to the surface 6 S . On the other hand, the larger density 
of s and p states at the surface and the rapid de- 
cay of the hybridization matrix elements with distance 
suggest that 6b is negligible, and the experiment for one 
impurity can also be explained if Vb = 0. A cal- 

culation of 6b has not been made and experimentally the 
situation is still unclear. The role of V& is not only crucial 
for a correct theory of the mirage experiment, but also 
for the general understanding of the interaction between 
metallic surfaces and adsorbates. Since actually 6b was 
introduced as a phenomenological parameter which takes 
into account the electrons lost in the scattering process, 
H[|Q] one expects that if 6b = 6 S , the interaction be- 
tween impurities is roughly a fourth of that for Vb = if 
the same total width 6b + 6 S is kept. 

We obtain the ground state of the Anderson model in 
a cluster which contains one or two impurities and the 
relevant conduction states inside a hard wall ellipse using 
the Lanczos method. These states are then mixed with 
bulk states using an embedding method. |l5| ] This em- 
bedding is essential to describe the low energy physics. 



]l4j The average separation between the relevant con- 
duction states d ~ 100 meV ||0] is much larger than 
the Kondo temperature Tr- ~ 5 meV. 0j Under these 
circumstances a Kondo peak at the Fermi level cf is ab- 
sent in the finite system and the experimental 
line shape for dl/dV cannot be reproduced (unless an 
artificial Lorentzian broadening is introduced). This is 
confirmed by our calculations. 

The Hamiltonian can be written as: 
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+ Y l V[ i p i [R i )4 tr c ia + H.c] + H'. (1) 

Here c] CT creates an electron on the j th conduction eigen- 
state of a hard wall elliptic corral with wave function 
<Pji r ) 01 an d d\ a is the corresponding operator for the 
impurity at site Ri. The hybridization of these states 
with bulk states of the same symmetry is described by 
H'. We assume that each of the impurity and conduction 
states mixes with a different continuum of bulk states: 
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H.c). (2) 



The bi a represent bulk states for which the unperturbed 
density is 0.05 states/eV, similar to the density of bulk 
s and p states. [0 Approximation (H) is justified by 
comparison of the non-interacting Green functions for 
hard wall corrals and more realistic boundary potentials. 
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The dressed matrix G describing the one-particle 
Green function is calculated by solving the Dyson equa- 
tion G = g 4- gH'G, where g is the corresponding ma- 
trix for H' = 0. 0] This equation is exact for U = 
or H' = 0, and in the general case represents an infinite 
sum of particular diagrams in perturbation theory in H' 
(the chain approximation [fig]). The <fij{r) are obtained 



as described elsewhere. ]1(J We choose the ellipse with 
eccentricity e = 1/2 and size such that the state j = 42 
lies at e F . @ The change in dl/dV (AdI/dV) after an 
impurity is placed inside the corral is determined by the 
conduction states which lie near tp and have a strong 
amplitude | (i?i)| at the impurity position. For Ri at 
one focus they are j = 32, 35, 42 and 51. We have 
also included j = 24 and 62, although this inclusion leads 
to negligible changes in the results. We took the impurity 
parameters Ed = — 1 eV and U = 3 eV. 0] We consider 
first the case Vb = and one impurity at the left focus 
{Ri = (-0.5a, 0)). The value V = 0.04 eV was chosen to 
lead to the observed width of AdI/dV. The remaining 
parameter t controls the amplitude of the mirage at the 
right focus. 

In Fig. 1(a) we represent the resulting impurity spec- 
tral density Pd{u) for two values of t. A clear Kondo 
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peak is obtained and for t > 0.3 eV its width is very 
weakly dependent on t. Instead, for t — > 0, the peak 
splits into two very narrow peaks out of tp. In contrast 
to the magnitude of the change in the conduction 

density Ap c (r, uj) at the empty focus (r = —Ri) is quite 
sensitive to t > 0.3 eV: as t increases, the width of the 
conduction states increases, the weight of the states 32, 
35 and 51 (odd under the reflection through the minor 
axis of the ellipse a) at ep increases, and the depression 
of p c (— Ri, uj) decreases as a consequence of the negative 
interference of these states with the even state 42. jl(J 
The differential conductance dl/dV at zero temperature 
is proportional to the density pf of the state (2(]] 

fa (r) = <Pj ( r ) C 3° + ■ ( 3 ) 
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q is related to Fano's interference parameter and repre- 
sents the effect of a direct tunneling from the tip to the 
impurity. Therefore, it is relevant only very near the im- 
purity. For q — 0, pf(r,u>) — p c (r,uj). In Fig. 1 (b) we 
represent the effect of adding the impurity on pj(±Ri, uj) 
(Apf ~ AdI/dV). At the impurity site Ri, Ap c (r,uj) is 
asymmetric and smaller at the right of the valley. This 
is a consequence of the asymmetry of the hybridization 
around ep (\ip^i(Ri)\ > | ^35 ) I ) • A symmetric line 
shape, as observed in the experimental AdI/dV is re- 
stored for q ~ 1. The effect of this q is consistent with 
the fact that on a clean surface, AdI/dV is larger at the 
right of the peak. Another nice fact is that the minimum 
of Apf for q = 1 lies at the experimental position 1 meV. 
At the right focus (r = — Ri) we obtain a similar valley, 
although slightly asymmetric and shifted to the left. In- 
creasing t from 0.4 to 0.5, the magnitude of this valley 
is strongly reduced (its minimum is shifted above -5/eV) 
but its shape and width is retained. At the impurity 
position there are no significant changes. 

The space dependence of Apf for q — is represented 
in Fig. 2. As in the experimental AdI/dV, the main 
features of | y?42 ("r) | 2 , attenuated at the right focus, are 
displayed. Thus, the theory reproduces the space and 
energy dependence of AdI/dV observed in the experi- 
ment. All results so far agree semiquantitatively with 
perturbative calculations. To see how the results 

change if 5b = S s is assumed, we have reduced V by a 
factor -\/2 • This should reduce 5 S by a factor 2. Increas- 
ing Vb from zero to 1.2 eV, the original width of pd is 
restored. The intensity is reduced by a factor ~ 2 (due 
to the strong energy dependence of S s ). Ap c turns out 
to be ~ 4 times smaller. The additional factor 2 can 
be understood from the fact that the change in conduc- 
tion electron Green function is proportional to V 2 Gd(ui), 
where Gd{uj) is the impurity Green function. ||To|| Except 
for these factors, the results are surprisingly similar to 
the previous ones. Some of them will be displayed in 
Fig. 4. 



We now turn to the case of two impurities, one at each 
focus, for Vb = 0. The spectral density for one of these 
impurities is represented in Fig. 3 (a). Comparison with 
the previous case (Fig. 1), shows that the peak around ep 
broadens (by a factor ~ 1.5), looses intensity and shifts 
to lower energies. In addition, another very narrow peak 
appears ~ 13 meV below ep. An analysis of the energy 
dependence of the density of the individual conduction 
states shows that the broad peak around cf is due to hy- 
bridization with even states (mainly 42), while the nar- 
row peak reflects the hybridization of the impurity states 
with odd states (mainly 51). The difference in dl/dV 
with respect to the empty corral is however, not so dif- 
ferent as in the previous case. This is due to the effect of 
the unperturbed Green functions of the conduction states 
and is also present in the one impurity case. (l(J Never- 
theless, a decrease in the amplitude and a broadening of 
the depression should be observed in AdI/dV and seems 
in qualitative agreement with recent experiments. Q The 
space dependence is similar to that for one impurity (Fig. 
2) but it is of course, symmetric under reflection through 
the minor axis a, and not attenuated at the right focus. 

Qualitatively, the shape of pd can be understood look- 
ing at the non- interacting case U — 0, Ed ~ ep. In 
this case, for one impurity, the Kondo peak is replaced 
by a Lorentzian near ep. For two impurities, a change 
of basis of the d orbitals to e a = {d\ a + duo) 
Oct = (filer — d2 a )/V2, separates the problem into those 
corresponding to even and odd states under a. The even 
state hybridizes mainly with conduction state 42, to form 
a resonance near ep, roughly twice wider than for one im- 
purity due to the larger effective hybridization. Instead, 
the odd state o a is displaced towards lower energies due 
to hybridization with state 51. The interactions should 
modify the quantitative details of this picture. However, 
we expect that it remains qualitatively valid, as suggested 
by the above mentioned energy distribution of the differ- 
ent conduction states. 

To gain insight into the nature of the ground state, we 
have also calculated spin-spin correlation functions for 
t = 0. A reliable method to include H' in these calcula- 
tions has not been developed yet. For one impurity we 
obtain (Si • S42) = —0.73, where Sj is the spin of the 
impurity i and Sj is the spin of the conduction state j. 
This value is close to the minimum possible one -3/4. For 
j 42, (Si • Sj) are very small, but this, and the large 
magnitude of (S, • S42), are affected to a certain degree 
by the neglect of H' in this calculation. The space de- 
pendence of (Si • s(r)), where s(r) the conduction spin 
at position r follows closely 1^42 (0| 2 - F° r two impuri- 
ties we find (Si • s 42 ) = -0.47 and (Si • S 2 ) = 0.21. In 
the limit of large U , one expects that the main features 
of the spin dynamics for Vb = are described by the 
Hamiltonian Hq = J(Si + S2) ■ S42, where J > is the 
Kondo coupling. The ground state of this Hamiltonian 
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is a doublet in which the impurity spins are correlated 
ferromagnetically between them ((Si • S2) = 1/4) and 
antiferromagnetically with state 42 ((Si • S42) = —1/2). 
These values are near to those we find. The effect of the 
hybridization of state 42 with bulk states can be mod- 
elled by a tight binding Hamiltonian in terms of Wilson's 
orbitals. Ho is the strong coupling fixed point of Wilson's 
renormalization group. An analysis of the stability of this 
fixed point using perturbation theory as in Ref. Q leads 
to the conclusion that the ground state is a doublet for 
Vb = 0. However, we expect that as soon as 14 7^ 0, the 
doublet is screened at a very low temperature. 

For the set of parameters corresponding to 5b = 5 S , 
Pd(u) is much more similar to the one impurity case, 
although a structure reminiscent of a splitting is also 
present near its maximum. In contrast to the case of 
Vb = 0, when a second impurity is added, the depres- 
sion in AdI/dV at one impurity site Rt increases and its 
width is roughly the same (see Fig. 4). Comparison with 
results when t is increased from 0.4 eV to t — 0.5 eV (not 
shown) suggests that when 8b = S s , AdI/dV at ± Ri for 
two impurities is roughly the sum of the results at Ri and 
— R4 for one impurity. This is what one would expect if 
the interaction is very small. 

Coming back to the case Vb — 0, we have also ver- 
ified that qualitatively similar features in AdI/dV are 
obtained at one focus, if one impurity is placed there 
and the second impurity is put at another extremum of 
<P42(?")> like (0.22a,0) (instead of placing it at the other fo- 
cus). In this case, the spectral densities at (0.22a,0) have 
some additional structure due to an important admix- 
ture of the state 41. In contrast, if both impurities 
are placed close to the same focus and near each other, a 
moderate hopping t ~ 0.15 eV or larger between them is 
sufficient to destroy the Kondo resonance. In particular 
AdI/dV becomes flat and featureless near cf- 

In summary, we have studied the spectral density for 
impurities inside a quantum corral, using a many-body 
approach which treats exactly the correlations in the im- 
purities and their hybridization with the relevant con- 
duction states at the surface, and treats approximately 
the hybridization with bulk states. We have been able 
to reproduce the main features of the mirage experiment 
for one impurity inside the corral. The experiment for 
one impurity cannot determine the relative importance 
of the direct hybridization of the impurity with bulk 
states, unless the tunneling matrix elements and other 
details are known accurately. Instead, for two impuri- 
ties inside the corral, the differential conductance is very 
sensitive to this hybridization. For the parameters of the 
experiment, the spins of both impurities are antiferro- 
magnetically coupled with the conduction electrons, and 
ferromagnetically correlated between them provided they 
are placed sufficiently far apart, so that the hopping be- 
tween them can be neglected. If this hopping is larger 
than 0.15 eV, there is a tendency to form a singlet state 



between both impurity spins and the Kondo resonance 
disappears. To our knowledge, this is the first theory 
which is able to describe the line shape of the differen- 
tial conductance when more than one Kondo impurity is 
inside the quantum corral. 
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Fig. 1: (a) Impurity spectral density as a function of 
energy for two values of t, (b) Change in the density 
of the mixed state f a (Eq.(|3|)) at the impurity site (left 
focus) for two values of q and at the other focus for t = 0.4 
eV. 

Fig 2: Contour plot of Ap c (r, u) for t = 0.4 eV and 
a; = 10 meV. 

Fig 3: (a) Impurity spectral density for one impurity at 
each focus and two values of t. (b) Change in the density 
of the mixed state f a after addition of both impurities 
(Eq.([|)), at one impurity site for two values of q. Param- 
eters are V = 0.04eV, V b = and t = 0.4 eV. 

Fig 4: Ap c (r, ui) as a function of u> for the case of one 
impurity at the left focus (full and dashed lines) or one 
impurity at each focus (dashed dot dot line). Parameters 
are V = 0.04eV/V2, V b = 1.2 eV and t = 0.4 eV. 
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